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A mathematical model is developed to analyze the growth/decay rate of surface per- 
turbations of an ultrathin metal film on an amorphous substrate {S1O2). The formulation 
combines the approach of Mullins [J. Appl. Phys. v30, 77, 1959] for bulk substrates, in which 
curvature-driven mass transport and surface deformation can occur by surface/ volume dif- 
fusion and evaporation-condensation processes, with that of Spencer et al. [Phys. Rev. 
Lett. v67, 26, 1991] to describe solid state transport in thin films under epitaxial strain. 
Modifications of Mullins model to account for thin film boundary conditions result in qual- 
itatively different dispersion relationships especially in the limit as kh <C 1 where k is the 
wavenumber of the perturbation and h a is the unperturbed film height. The model is applied 
to study the relative rate of solid state mass transport as compared to that of liquid phase 
dewetting in a thin film subjected to fast a thermal pulse. Specifically, we have recently 
shown that multiple cycles of nanosecond (ns) pulsed laser melting and resolidification of 
ultrathin metal films on amorphous substrates can lead to the formation of various types of 
spatially ordered nanostructures [Phys. Rev. B, v75, 235439 (2007)]. The pattern formation 
has been attributed to the dewetting of the thin film by a hydrodynamic instability. In 
such experiments the film is in the solid state during a substantial fraction of each thermal 
cycle. However, results of a linear stability analysis based on the aforementioned model 
suggest that solid state mass transport has negligible effect on morphology changes of the 
surface. Further, a qualitative analysis of the effect of thermoelastic stress, induced by the 
rapid temperature changes in the film-substrate bilayer, suggests that stress relaxation does 
not appreciably contribute to surface deformation. Hence, surface deformation caused by 
liquid phase instabilities is rapidly quenched-in during the cooling phase. This deformed 
state is further evolved by subsequent laser pulses. These results have implications to devel- 
oping accurate computer simulations of thin film dewetting by energetic beams aimed at the 
manufacturing of optically active nanoscale materials for applications including information 
processing, optical devices and solar energy harvesting. 
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I. INTRODUCTION 

Ion and photon beam irradiation are widely used to generate nanoscale surface features for a va- 
riety of applications: optoelectronic, plasmonic, magnetic, bio- and chemical sensing, and catalytic 
devices [3,13, 0,0] • Understanding interfacial transport mechanisms under fast thermal processing 
or energetic beam irradiation is essential to accurately predicting the spatio-temporal dynamics 
of film evolution that imparts the desired nanoscale surface features. The contributions from var- 
ious solid state transport processes that cause surface changes such as flattening or grooving in 

nnQ 

bulk materials were studied by Mullins [a, la, lD| • In his formulation, smoothening of bulk surfaces 
is driven by curvature effects that manifest in solid state surface diffusion, volume diffusion, or 

n 

evaporation-condensation mechanisms [7||. Mullins' theory has been widely accepted and used ex- 
tensively in the analyses of surface processes Q, 0, Q, 11, Q, 3|. However, Mullins' theoretical 
framework is inadequate to describe solid state mass transport in ultrathin films in which one has 
to account for the effects of other driving forces on pattern formation, e.g. dispersion forces and 
epitaxial stress. For instance, it is well-established that epitaxial stress has a significant effect on 
surface diffusion. There is extensive literature on pattern formation in thin films in the presence of 



epitaxial stresses [1J, LL 
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1 71 ] - In contrast, the literature on the evolution of ultrathin films on 
amorphous substrates, a scenario in which epitaxial stress effects are negligible, is sparse. In this 
article, we develop a generally applicable mathematical model for describing such processes and 
apply the model to study solid state transport in thin films exposed to time periodic laser thermal 
stimulus. 

Recently, we have observed that thermal processing of ultrathin metal films by multiple instances 
of a 9 ns pulsed laser (with each pulse separated in time by 20 ms) results in a variety of spatially 
ordered patterns Q, Q, O ne of the important findings, relating to the processing conditions, 
was that in order for pattern formation to occur, it was necessary to irradiate the films at a laser 
fluency at or above a critical value, which was well-correlated with the threshold energy required to 



melt the films [211 ] . Consequently, mass transport was considered to occur primarily in the liquid 
phase and pattern formation was attributed to a thin film hydrodynamic (TFH) dewetting instabil- 
ity, which arises in spinodally unstable films when attractive long-range intermolecular dispersion 
forces amplify surface perturbations [^J. Thermal modeling of the film/substrate bilayer showed 
that during one laser pulse, the film experienced a nearly uniform, but transient temperature field 
for a time scale of 0(100 ns), of which the time to reach the melt temperature during the heating 
phase was ~ Ins, the liquid lifetime, T£ was typically between 1 — 10ns and cooling time was 
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~ 100 ns [21[. Therefore, approximately 90% of the thermal cycle following each laser pulse was 
spent in the cooling phase. This is important in the context of Mullins' analysis, which suggests that 
solid state diffusion and evaporation-condensation processes could influence surface deformation. 
Hence, it is necessary to analyze both liquid and solid phase mass transport processes to develop 
a comprehensive understanding of the physics of interface evolution in thin films processed under 
pulsed laser irradiation. 

In this work, we have suitably modified Mullins' original formulation j?J] to describe the evolution 
of infinitesimally small surface perturbation of ultrathin solid films on amorphous substrates. In 
particular, we have derived new results for volume diffusion transport and have employed the 
approach adopted by Spencer et al. [14j] originally for describing diffusion in epitaxially stressed 
films, in order to account for dispersion forces that are influential only in thin films. For pulsed laser 



experiments described in 



or cusp 
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2 If ] , we determined that there is negligible contribution to surface 
deformation from solid state material redistribution by evaluating solid state mass transport rates 
via diffusion and evaporation-condensation processes. These results suggest that film deformation 
occurring under each melting cycle in the liquid phase is essentially quenched-in during rapid 
resolidification. Subsequent laser pulses cause progressively more interface 
film ruptures. This mechanism is consistent with experimental observations 
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ii. SOLID STATE MASS TRANSPORT 

Local curvature differences provide the driving force for mass rearrangement through vari- 
ous mechanisms in the solid state including surface diffusion, volume diffusion, and evaporation- 
condensation. In a series of papers, Mullins outlined a theory to predict the mass transfer rates 
by the aforementioned processes for surface deformation of bulk materials . For thin films, 

additional driving forces need to be accounted for, such as the free energy contributions from dis- 
persion (van der Waals) forces and thermal strain. For films supported by amorphous substrates, 
epitaxial strain is irrelevant. However, mismatch in thermal expansion coefficients of the film and 
the substrate can result in thermoelastic strain. In this work, we first develop the theory by incor- 
porating dispersion forces, in the absence of thermal stresses. The effect of thermal stresses is only 
discussed qualitatively in Se c. [Vl 

Following Spencer et al. [LJ], the evolution of the film height h is described by: 

dh 

— = B SD V 2 ( 1 K + e) + V (1) 
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where B$d = DsQ 2 v/ksT, 7 is the surface tension, k is the local curvature, D$ is the surface 
diffusion coefficient, O is the atomic volume, v is the number of atoms per unit area, fee is Boltz- 
mann's constant, T is the absolute temperature, e is a free energy per volume term, and V is an 
additional velocity term such as a material deposition term jj], Q|. In the limit that the height of 
any surface fluctuations are small, the curvature can be approximated as k ~ V 2 /i. It has alr eady 



been established that dispersion forces play a 



especially upon laser-induced melting 
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e in the evolution of a thin metal film [221 ] . 
23(| . By including the effect of the disjoining 



pressure to the free energy, the height evolution equation becomes: 

% = B SD V* b V> h+ ^ ]+ V (2) 

where A is the Hamaker coefficient and the free energy per volume contribution due to dispersion 
forces, g^p-, is strictly valid for a flat film. 

The effects of evaporation-condensation and volume diffusion processes can be included into the 
height evolution equation as additional velocity terms, in eq. [TJ By relating local pressure differences 
between the flat and perturbed states to the local curvature though the Gibbs-Thompson formula, 
Mullins has shown that the evolution of a rough surface, due to evaporation-condensation is given 
by: 

dh 

Vec = (-^)ec = B EC V 2 h (3) 

where Bec = PnyQ 2 /(2nm)2 (k B T)% , p a is the equilibrium vapor pressure and m is the mass of the 
metal atom [j], Q|. 

As in the case of evaporation-condensation, volume diffusion can be considered and included in 
the height evolution equation as another V term. In the derivation below, x and z denote the lateral 
and normal coordinates respectively such that z = and h correspond to the film-substrate interface 
and local film height respectively. The substrate is much thicker than the film and hence considered 
as a semi-infinite domain. The decay rate due to volume diffusion can be related to the local density, 
p as H = DyQ(^) z= h, where Dy is the coefficient of self-diffusion and z denotes the coordinate 
normal to the substrate plane F° r a bulk surface (—00 < z < 0), assuming quasi steady-state 

diffusion (V 2 /?(x, z) = 0) and solving the Laplace equation with boundary conditions p(x, z) = p a 
at z — > —00 and p{x,z) = p + p -^pk 2 esin(kx) at z = 0, where e and k are the amplitude 
and wavevector of the surface perturbation h' respectively, Mullins derived = Byoik^h', where 
By i = D ftjr ■ However, for a film of finite thickness the boundary conditions need to be modified, 
specifically: p(x,z) = p Q at z = and p(x,z) = p a + p -^^k 2 tsin{kx) at z = h. Assuming quasi- 
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steady state diffusion and h = h Q + h! where h' = ee ast sin(kx) the expression for volume diffusion 
becomes <C 1): 

Vv ° l = ^ = ~ Bvolk 1 - exp(-2kh o) h (4) 
where h Q is the average film thickness. 

h 2 

Note that the quasi-steady state (QSS) assumption is only valid if t -p 2 -, and even for sub- 
nanometer thick metal films, the time to reach QSS, t S> 200ns. However, since (^) 2= /j (the 
density gradient at the surface of the film) will have its largest magnitude for short times and 
decrease to its steady-state value, the decay rate represented by a quasi steady-state approximation 
is a limiting case. For systems with a timescale t < jy- the decay rate will be faster than that 
predicted above. For macroscopic systems {kh » 1), eq. |4] reduces to the Mullins equation for 
volume diffusion, dh' /dt = —By^k^h' [7]. It is interesting to note that in the opposite limit, i.e., 
for thin films with a relatively large perturbation wavelength {kh <C 1) and after linearizing the 
exponential about h Q , the expression for decay due to volume diffusion simplifies to: 

dh' -B Vol k 2 

-dt=—h^ h - (5) 

This limit applicable for thin films is qualitatively different from that for a bulk material, with the 
rate of decay varying as k 2 as compared to k 3 for the latter. Further, the thickness dependence, i.e., 
the decay rate varying as is absent in the case of bulk materials. 

By combining eqs. [21 [3] and HI the evolution equation for height perturbations h'(<^ h Q ) can be 
expressed as: 

Letting h! = ee ast sin(kx), the growth/decay rate, o~ 5 of surface perturbations in the solid state is 
given by: 



r r ;4 A 7„2i o ,„3 1 + exp(-2kh Q ) 2 
as = -Bsohk - ^k] - B Vol k - 1 _ exp( _ 2kho) ~ B EC k (7) 



III. HYDRODYNAMICS OF THIN LIQUID METAL FILMS 



In this section we summarize the results for the growth/decay rate of thin liquid films with the 
purpose of comparing the mass transport rates in the liquid state to that in the solid state. For 
a detailed description of liquid state mass transport and dewetting in thin films, see [2a, |2a, |27j, 
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28|, |29|, |30| . When the Navier-Stokes equations, describing mass transport in thin liquid films, are 
simplified via the lubrication approximation, a film-thickness dependent relationship between the 
material flow rate and pressure gradient can be derived 



25 



261 ] . The pressure is a combination of the 



Laplace pressure from capillarity and a disjoining pressure that arises from long-range dispersion or 



van der Waals forces 



25], 



261 ] . When the disjoining pressure, which tends to amplify perturbations, 
exceeds the stabilizing effect of the surface tension, the film is in an unstable state and perturbations 
will grow, eventually causing the film to rupture. The thin film hydrodynamic equation describing 
this behavior is given by the following: 

^ = __L V ■ ( 7 h 3 VV 2 h + -Ay/*) (8) 
dt 3rj w 2irh 1 v ' 

where r\ is the dynamic viscosity. Once again letting h = h Q + ee aLt sin(kx), the growth/decay rate 
CTiin the liquid phase is given by: 



yielding the largest unstable wavevector, k max = * / - — and a fastest growing perturbation wave- 

V ^l^o 

length, A* = J^hl 



IV. EVOLUTION OF SURFACE PERTURBATIONS IN SOLID AND LIQUID STATES 



As mentioned previously (Sec. IIH|) . analytical models and numerical simulations have shown 
that for ns pulsed laser irradiation of thin metal films, the cooling time t coo i of the film is approx- 



imately 10 times as large as the time that the film is in the molten state, i.e., t coo \ ~ 100ns [21]. 
At the beginning of the cooling phase, the film is at its melting temperature. The high temper- 
atures suggest that solid state material transport mechanisms could be active. To estimate the 
relative magnitude of the effect that solid state mass transport could have on deformation of the 
film during the cooling period, we employed eqs. [7] and M For the values of D$ and Dy, we used 
the models of Flynn for surface and volume diffusion in which D$ = 5 x 10 _4 exp(— 6Tm /T) 



and Dy = 0.3exp(— VJTm/T) [3l|]. In order to gain quantitative insight on the effect of mass 
transport in each phase, we have compared the growth/decay rate of a perturbation multiplied by 
the timescale during which it is applicable. Since diffusion and evaporation processes are dispro- 
portionally influenced by the time spent at the peak temperature, we have used an overestimate of 
the time at the peak temperature of Tg = 50 ns in order to find an upper limit of the contribution 
from the solid state transport mechanisms. In Fig. [T^a-b) we have plotted the magnitude of as *t$ 
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and ai, * rjj as a function of perturbation wavevector k for material redistribution in both phases, 
where r| = 10 ns is the typical liquid lifetime. Due to the thickness dependency of several of these 
mechanisms, we have compared transport mechanisms in a lnm and a 10 nm Co film. 

In the case of a 1 nm film, the stability of the solid film is strongly dependent on the choice of 
the Hamaker coefficient A, which dictates the magnitude of the disjoining pressure. When using the 
experimentally determined value of A = 1.4 x 10 -18 J {^(J in eq. [3, the film is unstable and there 
exists a characteristic wavelength that will grow the fastest. The dispersion curve for growth/decay 
rates of perturbations in the solid state has been shifted to larger wavelengths, as compared to that 
for the TFH instability, as shown in Figure (Ha). In this case the cut-off and the fastest growing 
wavelengths have been increased due to the stabilizing effects of volume diffusion and evaporation- 
condensation. In the absence of of these effects k max and A* are the same for both the solid and 
liquid states. The importance of the magnitude of the dispersion term on mass transport can be 



3J. When a 



seen if one uses the lower limit of the Hamaker coefficient for metals (~ 10~ 19 J) 
theoretically estimated A ~ 3 x 10~ 19 J was employed, the disjoining pressure was not large enough 
to destabilize the solid film; but rather, the film was seen to be stable against perturbations of all 
wavelengths. For a 10 rem film, the dispersion forces with A = 1.4 x 10~ 18 J are not large enough 
to instigate a dewetting instability in the solid film, as shown if Figure [0(b) . For the length scales 
stemming from a TFH instability, | &L |^>| as |, thereby confirming that hydrodynamic flow is the 
more dominant mass transport mechanism during ns pulsed laser melting. We obtained similar 
results for a variety of metals suggesting that this is a general result, applicable to most thin metal 
films. Most metals investigated have similar diffusion rates, but the evaporation rates vary greatly. 
Figs. [2la-b) show plots of dispersion curves for two extreme scenarios, one metal (Al) with a low 
evaporation rate and the other (Cr) with a large rate. In both cases, films with h Q > 2nm are seen 
to be stable to small perturbation of all wavelengths. 

There are several implications manifested in these results. First, the timescales of solid state 
process are much larger than the timescales associated with hydrodynamic flow. Therefore, the 
quenched-in structures from the pulsed laser melting and resolidification are likely unchanged dur- 
ing the cooling period. This conclusion is corroborated by experimental observations of pattern 



formation in ultrathin metal films by using multip 



ultraviolet (UV) laser irradiation: see [1 



21 



e pulses of (266 nm wavelength, 9 res pulsed) 



331 ] for details. A variety of metal thin films, 



including Co, Fe, Cu, Ag, Ni and Pt, were deposited under ultrahigh vacuum (~ 10~ 8 Torr) by 
e-beam evaporation or pulsed laser deposition onto commercially obtained SiC^/Si substrates with 
the thermally grown oxide layer of 400 nm thickness. Film thicknesses ranged from 1 — 10 nm. The 
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initial films had an average surface roughness < 0.2 nm as evaluated by atomic force microscopy. 
Following deposition, the films were irradiated at various laser energies above a threshold energy 
density E t h such that appreciable surface deformation was observed for various irradiation times, 
as measured by the number of pulses n ranging from 10 < n < 10, 500. Following such irradi- 
ation, the surface morphology and concentration of metal was evaluated by imaging and energy 
dispersive X-ray spectrometry in a scanning electron microscope, respectively. It was found that all 
the metals possessed a thickness-dependent threshold energy density E t h above which substantial 
surface changes were observed. However, for energies below this threshold, no surface roughening 
was visible even after the longest irradiation times of 10,500 pulses. This energy threshold has been 



211 ] . By irradiating films with a laser energy 



attributed to the energy required to melt the films 
density above the melt threshold, we performed detailed studies of the dynamics of the pattern 
evolution for Co films as a function n and found that the dewetting patterns showed characteris- 
tic film-thickness dependent length scales and spatial order consistent with the predictions of eq. 
[9] Hence, despite the much larger time spent in the solid phase (i.e. cooling period), the overall 
contribution from solid state transport was inferred to be negligible. 

V. THERMAL STRESS 

Since we focus on the dynamics of perturbed thin films on amorphous substrates, the prevalent 
stresses that are present in the film are thermally generated. Due to the mismatch in the thermal 
expansion coefficients between the film and the substrate, the film is under compressive stress while 
it is heated and a tensile stress when it is cooling. This stress is given by the following expression: 

*/ = ( °' = a/) f T (10) 

1 ht\ s + ±=g*l 

n nh s E s t h a E a \ 

where a is the thermal expansion coefficient, v is Poisson's ratio, E is Young's modulus, h is the 
thickness of the layer, AT is the temperature change, and the subscripts / and s denote the film 
and substrate quantities, respectively [34]]. In the limiting case where j^f- S> \ f _^J , the preceding 
expression reduces to: 

a f = ^"^ AT. (11) 

Stresses related to thermal expansion for most thin metal films on a silica substrate are on the order 
of a few GPa. However, they are still lower than the critical value required to initiate a buckling 
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instability in the film, which is given by the following equation: 

Ef i E s .2 

^critical = Hz — ) 3 ( , _ 2 ) 3 (12) 

f V& 

where 4> = 0.52 is a constant. Typical critical stress values range from ~ 40 — 60 GPa. Stress 



can also result from an anisotropic temperature distribution in the film [35J]. This could arise in 
pulsed laser processing. As described above, the thickness variations caused by the TFH instability 
are quenched in during the rapid resolidification phase. Hence, in the solidified deformed film, 
a non-uniform temperature profile could prevail for times comparable to the diffusion time scale 



211 ] . Within the context of linear theory of elasticity, the thermoelastic stress will also decay at 



351 ] . Hence a characteristic time scale 



the same rate as the the temperature profile in the film 
of temperature/stress diffusion can be estimated by solving the transient heat equation for a 2- 
dimensional film: 

„ 8T 8 2 T d 2 T , 
P ' Cp J-dt = K 'W + « f d* (13) 

where pf, Cpj, and Kf are the density, heat capacity and thermal conductivity of the film, respec- 
tively. Assuming an initial sinusoidal temperature profile in the lateral (x) direction and a uniform 
temperature distribution throughout the thickness [z direction) of the film, one finds the following 
relationship for the time dependent component of the temperature, T'(t), in which the temperature, 
T(x,t)=T (x)T'(t): 



T'(t) = exp[-a f kH - ^f Vi] (14) 

Pf^PJ^o 

ctf is the thermal diffusivity of the film, k is the wavevector of the perturbation, c is a constant 
between and 1 and p s , Cp )S , and k s are the density, heat capacity and thermal conductivity of the 
substrate, respectively. Since the length scale of the periodic thickness variations is being imprinted 
on the film in the liquid state, one can take the expected value of k for a corresponding h from 
TFH theory. Using these values (c = 0.2), the timescale during which this stress could influence 
solid state mass transport is on the order of a picosecond for a perturbation wavelength, A ~ 30 nm 
and on the order of a nanosecond for A ~ 1 pm. These timescales are too short to appreciably effect 
any solid state process in modifying the surface of the film. 



VI. CONCLUSION 



We have developed a mathematical model to describe thin film dynamics on amorphous surfaces 
by incorporating the effects of surf ace/ volume diffusion, evaporation-condensation and dispersion 
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forces. New results are presented for volume diffusion in thin films. Specifically, for kho< < 1, the 
decay rate of surface perturbations by this mechanisms is proportional to k 2 /ho as compared to 
k 3 for bulk materials. It is shown that depending on the magnitude of the Hamaker coefficient, 
ultrathin solid films with thickness < 2 nm could be prone to pattern forming instabilities while 
films with thickness > 2 nm are predicted to be stable for a variety of metals. The model was 
applied to analyze the role of solid state mass transport mechanisms on the dynamics of thin films 
processed by pulsed laser irradiation. The upper limits for the time constants associated with the 
decay of surface perturbations in the solid state was found to be much smaller compared to the 
time interval between subsequent pulses. Further, a qualitative analysis of thermoelastic stress 
dissipation in the film-substrate bilayer system showed that the time scale of stress relaxation is 
much smaller than the pulse duration. These findings suggest that contributions to the changes 
in surface deformation from solid state processes are minimal and the pattern formation occurs 
through a hydrodynamic instability that is followed by laser melting of the film, an observation 
corroborated by experimental observations. 
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Figure 1: Plots of the dispersion relations for (a) 1 nm Co film (using A = 1.4 x 10~ 18 J), with the scale 
factor, X\ — 2000 and (b) 10 nm Co film where scale factor X w — 100. In both plots Tg = 50 ns and 
t p l = IQns. 
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Figure 2: Plots of the dispersion relations for (a) Aluminum and (b) Chromium solid films with thicknesses 
of 1 nm and 2nm, using and estimated Hamaker coefficient A = 5 x 10~ 19 J. Both types of films are stable 
to small perturbations of all wavelengths for h >^ 2nm. 



